Matter Wave Scattering from Ultracold Atoms in an Optical Lattice 
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We study matter wave scattering from an ultracold, many body atomic system trapped in an 
optical lattice. We determine the angular cross section that a matter wave probe sees and show 
that it is strongly affected by the many body phase, superfluid or Mott insulator, of the target 
lattice. We determine these cross sections analytically in the first Born approximation, and we 
examine the variation at intermediate points in the phase transition by numerically diagonalizing 
the Bose Hubbard Hamiltonian for a small lattice. We show that matter wave scattering offers a 
convenient method for non-destructively probing the quantum many body phase transition of atoms 
in an optical lattice. 



Ultracold atoms in an optical lattice form an optical 
crystal that permits the clean implementation of funda- 
mental models of condensed matter physics to a degree 
unheard of in other systems. The physical parameters 
that characterize this system, namely the depth of the 
lattice and the interactions of the atoms in the lattice, are 
experimentally tunable so that it is possible to probe the 
atoms' behavior controllably from non-interacting single 
particle mechanics to strongly interacting, many body 
physics [1,2.3. Of particular interest are the many body 
theories of quantum phase transitions, whose predictions 
can be observed in such a system. 

The superfluid to Mott insulator phase transition has 
been probed experimentally by studying the interference 
of clouds of atoms expanding freely out of an optical lat- 
tice [1]. This method provides information about the cor- 
relation properties of the ground state of the cold atom 
system, but does not depend on the excitation spectrum 
of the lattice, which is essential to study the superfluid 
fraction [5]. Bragg spectroscopy has been combined with 
this method in order to examine the excitation spectrum 
[BJ; nevertheless, it requires the destruction of the sample 
under examination. More recently, it has been proposed 
theoretically that the on-site number statistics can be 
probed by observing light scattered into an optical cav- 
ity [7] . Very recently, light scattering from optical lattice 
systems has been analyzed to determine the effect of the 
many body phase on the scattering cross section of a 
photon due to interactions with a lattice [8]. Very lit- 
tle attention has been paid to the interaction of matter 
waves with these optical crystals. Theoretical work has 
been done emphasizing the impact of disordered atoms 
in a periodic potential on scattering of photons and par- 
ticles [9]. We are interested, however, in the impact of 
the many body phase on matter wave scattering. 

The scattering of matter waves from a lattice system 
provides a very useful technique to probe the many body 
phase transition both because it does not require the de- 
struction of the lattice under examination and because 
it depends strongly on the excitation spectrum of the 



target. This is critical to probe superfluidity and the 
dynamics of the lattice, beyond analysis of ground state 
properties. In addition, the simple form of the interac- 
tion between a slow-moving probe atom and the atoms 
in the lattice emphasizes structure that depends on the 
many body properties of the lattice. As we will show, 
the differential cross section is a highly suitable quan- 
tity for the distinction between the Mott insulating and 
superfluid phases. 

Our objective has been to study the scattering patterns 
of matter waves due to passage through ultracold atoms 
trapped in a uniform optical lattice. The scattering tar- 
get is well modeled by the Bose Hubbard Hamiltonian, 

H B h = -J ^2 + ^ U< y] n R (h R - 1) (1) 
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The probe is a free particle of mass, m, with Hamiltonian, 
Hp = p 2 /2m, which does not interact with the confining 
light of the lattice. We require low-energy probes that 
avoid intcrband excitations of atoms in the lattice. In this 
case, s-wave scattering is dominant, and we may treat the 
interaction between the probe and each lattice atom as 
a pseudopotential with scattering length, a s . The total 
interaction potential is V = Y^j ^r a sH^~ r?) EH- The 
operators r and Tj give the positions of the probe and 
the j th lattice atom, respectively. The full Hamiltonian 
for the scattering interaction is H = Hp + Hbh + V. 

The states of the lattice before and after scattering 
of a probe atom are many body states of the N atoms 
in the lattice. The particular ground state in the lat- 
tice depends on the relative sizes of the the interaction 
strength, U, and the tunneling matrix clement, J, ap- 
pearing in the Bose Hubbard model. For weak repulsion 
between the atoms in the lattice, the atoms will delocal- 
ize and the superfluid fraction will increase to one as the 
interaction strength goes to zero [5]. As the repulsion be- 
tween the atoms in the lattice becomes large compared to 
the tunneling matrix element, the atoms will localize, the 
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superfluid fraction will go to zero, and a gap will open in 
the excitation spectrum, giving rise to the Mott insula- 
tor state. It is possible to alter the interaction strength 
between the atoms in the lattice by adjusting the depth 
of the lattice, or by manipulating the scattering length of 
lattice atom collisions through a Feshbach resonance. It 
is best for the purpose of probing the many body phase 
of the lattice to retain a constant lattice depth so that 
the scattering patterns are not trivially affected by the 
changing density profile associated with a changing lat- 
tice potential. 

The incident probe wavefunction is in the plane wave 
state, |feo), with energy, ti 2 k 2 / (2m) . The target is ini- 
tially in the ground state, \n ), of the Bose Hubbard 
Hamiltonian, with energy, E no . We consider the cross 
section for a transfer of momentum k = ko — k from the 
probe to the lattice. This is associated with a change 
of the state of the atoms in the lattice to a poten- 
tially excited state, |n). The energy gained by the lat- 
tice atoms is related to the energy lost by the probe by 
E n — E no = j— (&Q — k 2 y In the first Born approxi- 
mation, the scattering cross section separates into two 
factors. One of which is determined by the binary inter- 
action between the probe and each target. The other is 
determined solely by the structure of the target [TTJ . This 
is true for general choices of the interaction potential. In 
the case of the pseudopotential, the scattered wave from 
an individual target is a structureless spherical wave, and 
the total angular cross section is given by 



the ground state of the pure superfluid target is 
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The cross section vanishes when the quantity under the 
square root becomes less than zero, due to energy con- 
servation. The momentum transfer, k, depends on the 
index, n, of the final state through the dependence on 
the final wavenumber, k, of the probe. A natural choice 
of units for the cross section is the square of the scat- 
tering length, a s . The quantity under the square root 
diminishes the contribution of scattering into final target 
states at progressively higher energy. The matrix element 
of the momentum boost gives the transition amplitude 
between the ground and excited states of the target due 
to a transfer of momentum, n, from the probe. When 
the interaction strength between the lattice atoms, U, 
is very large or small compared to the tunneling matrix 
element, J, it is possible to evaluate the cross section 
in ([2]) analytically. At intermediate values of U/J, we 
will numerically diagonalize the Hamiltonian to evaluate 
the cross section. 

For very weak interactions, we may treat the target as a 
condensate in the lowest energy Bloch wave, ipo(r), of the 
lowest band of the lattice. Designating the operator that 
creates particles in this state by = Jd 3 rijjo(r)^(r), 
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TV is the number of atoms in the lattice, Nl is the number 
of lattice sites. In order to calculate the matrix element 
appearing in ([2]) for momentum transfer, n, to the target, 
it is useful to express the many body operator in terms 
of the density, n(r) — (r)-0(r), using Ylj=i e%K rj = 
j d 3 r e ZKr n(r). 

A convenient basis for the target in the non-interacting 
case is a designation of the number of atoms in each mode 
of the lowest band of the lattice. Matrix elements of the 
form, (n\ n(r) \n ), are non-zero if at most one atom in 
the condensate has been excited out of the lowest energy 
Bloch wave. We only need to consider final states of the 
form, |n) = tp q where |£) represents N — 1 atoms in 

the ground state and ip q creates particles in the q ^ 
mode. The matrix element then takes the form 



(n\ n(r) \n Q ) 



(4) 



The energy difference of the target state, when only a 
single atom is excited, reduces to the single particle en- 
ergy difference between the two Bloch waves, E n — E nn — 
e(q) — e(0), where e(q) is the energy of a single atom in 
the specified Bloch mode of the lattice, q, corresponding 
to the excitation. The \n) — \n ) case must be handled 
separately. The result we obtain for the diagonal matrix 
clement is 



(n \ h(r) \n ) 



N 9 
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(5) 



We may now write down an expression for the cross 
section in ^ in the case of a superfluid ground state us- 
ing the above results for the matrix element of the density 
operator, 
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The q = component of the cross section gives Bragg 
peaks due to the coherent overlap of elastically scattered 
waves from each individual lattice site. The scale of the 
elastic scattering that gives rise to the Bragg peaks can 
be determined by considering scattering in the forward 
direction (n = 0). There we find that the central Bragg 
peak has a height ^ — N 2 . 

The sum in the second term is over the modes of the 
lowest band of the lattice. This term is exactly N times 
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the single target cross section. This includes the elastic 
single target channel, q = 0. We can estimate the scale 
of this term by considering the case in which the probe 
energy significantly exceeds the bandwidth of the low- 
est band, but is insufficient to excite atoms into higher 
bands. The energy of the probe and depth of the lat- 
tice are conveniently specified in units of the recoil en- 
ergy, E r — h 2 k\ I (2my) , for photons with wavenumbcr 
and lattice atoms of mass, tut- This condition is read- 
ily achieved for a typical lattice depth of V a = X5E r , in 
which the width of the lowest band is 0.03E 1 ,- and the 
band gap between the first and second bands is 6.28E r . 
When we approximate the final wavenumber of the probe 
to be equal to the initial wavenumber, the second term 
in Q becomes Na 2 .. 




FIG. 1: The analytic results for the superfluid (light gray) 
and Mott insulator (dark gray) differential cross sections for 
a sample lattice in ID with 15 atoms and 15 sites. These 
are shown on a log scale in order to draw attention to the 
absence of inelastic scattering from the Mott insulator. The 
cross section is for a probe with energy 6E r and a lattice 
depth of Vo/E r = 15. 

There are two major features to the scattering from 
a superfluid: narrow elastic Bragg peaks that scale as 
iV 2 and a superimposed inelastic background that scales 
as N. Fig. [T] which shows the differential cross sec- 
tions for the superfluid (light gray) and the Mott in- 
sulator (dark gray), illustrates this behavior for a one- 
dimensional lattice arranged perpendicular to the inci- 
dent probe wavevector. We consider the angle of devia- 
tion in the plane of the lattice. As the number of lattice 
sites increases, the width of the elastic peaks will be- 
come increasingly narrow. Away from the sharp Bragg 
peaks, the inelastic background is readily identifiable. 
This background emerges due to the availability of ex- 
cited state modes to the condensate. The scattering 
behavior is qualitatively different when the interaction 
strength between atoms in the lattice becomes very large. 



As the atoms in the lattice repel each other more 
strongly (U/J — > oo), the superfluid fraction decreases, 
and the atoms become localized within individual wells 
of the lattice. For a sufficiently strong interaction, the 
Mott insulator state forms, and the ground state of the 
target can be represented by the number of atoms at 



each lattice site, 



A uniform lat- 



i = n Rl ,...,n Rh 

tice will have n = N/N^ atoms per site for integer n. We 
must determine the matrix element (n\ h(r) \n Q ) for this 
ground state. This is most easily done by expanding the 
density operator in a Wannier basis, using the Wannicr 
function for the lowest band of the lattice, w(r). Then 
H r ) = Y,r u r 2 w*(r-Ri)w(r-R2)a Ri aR <i . Each term 
in the sum gives the contribution of the process in which 
a single target atom is scattered from one site to another 
by the probe. The R± = R2 term is the elastic channel 
in which the state of the target is unchanged, |n) = \n ). 
Inelastic scattering corresponds to R\ 7^ i?2. The ma- 
trix element of the final state of the target, in which 
one atom has been displaced from R' to R, is given 
by (n\ h(r) \n ) = ^(h+l)nw*(r - R)w(r - R). The 
energy cost associated with displacing one atom from 
the uniform ground state is the interaction strength, U, 
which is also the size of the Mott insulator gap. These 
results permit us to write the explicit expression for the 
scattering cross section of the Mott insulator target, 
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The sum over lattice sites in the first term takes a max- 
imum value when the momentum transferred from the 
probe is a reciprocal lattice vector, with k ■ R an integer 
multiple of 2n. As the number of lattice sites increases, 
the elastic Bragg peaks become increasingly sharp. In 
addition, there is an approximately Gaussian envelope 
due to the Fourier transform of the Wannier function. 

We can examine the k = case of forward scattering, 
as we did for the superfluid cross section. We see that 
the central peak has a height = iV 2 . The elas- 

tic scattering from the Mott insulator overlaps strongly 
with the elastic scattering given by the superfluid; how- 
ever, inelastic scattering from the Mott insulator phase 
is strongly suppressed (see Fig. [TJ. In particular, if the 
incident energy of the probe is less than the Mott insula- 
tor gap, the inelastic scattering vanishes completely. For 
probe energies exceeding the gap, the inelastic scattering 
is also negligible under the tight binding approximation. 
In that case, the integral in the inelastic part of the cross 
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section in Q is negligible when R ^ R', so that we 
expect only elastic scattering from the Mott insulator. 

We can estimate the scale of the Mott insulator's 
inelastic background more precisely by using the har- 
monic approximation of the Wannier function, in which 
we substitute the ground state of the harmonic ap- 
proximation to the bottom of an individual well in the 
lattice. Near to the central peak, the inelastic back- 
ground decays exponentially with the lattice strength, 

as exp (^—\ \p^-{j — 2 ) i where j and I (j ^ I) are po- 
sitions of lattice sites in units of lattice spacings. For 
lattice depths in the typical range V Q < 30E r , the inelas- 
tic contribution is strongly suppressed even for incident 
probe energies that exceed the gap energy. This contrasts 
markedly with the superfluid cross section, which carries 
a prominent inelastic background that scales as the num- 
ber of atoms in the lattice. In the regions between the 
Bragg peaks, this background serves as an unambiguous 
indicator of the many body phase of the lattice. 

We have also examined the disappearance of the su- 
perfluid inelastic background as the interaction strength 
between the lattice atoms is increased. This required that 
we calculate scattering cross sections for arbitrary values 
of the parameter U/J. At intermediate values, this re- 
quires knowledge of the spectrum of the Bose Hubbard 
Hamiltonian. We calculated the matrix elements of the 
density operator by exactly diagonalizing the Bose Hub- 
bard Hamiltonian for small lattices. Using these results, 
we have shown the angular cross section given in (§ for 
several values of U/J in Fig. [2] At U/J = the nu- 
merical result coincides with the analytic result we pre- 
sented. The inelastic background that scales as the num- 
ber of atoms is present. This inelastic background decays 
to zero as the interaction strength is increased, and the 
cross section converges on the analytic result we gave for 
the Mott insulator. We note that the amplitude of the 
inelastic background has decayed by more than half at 
U/J=4. For U/J=8, the background is largely gone. 
This coincides with the range over which the superfluid 
fraction vanishes [5]. 

Our analytic results for the scattering cross section 
when the interaction strength dominates the tunneling 
matrix clement (U J) and vice versa show that the 
many body phase in the lattice strongly affects the scat- 
tering cross section of a low-energy matter wave probe. 
The periodic nature of the target gives rise in both many 
body phases to coherent Bragg peaks whose height scales 
as the square of the number of atoms in the lattice. In 
addition, an inelastic background determined by the ex- 
citation spectrum of the target serves as an indicator of 
the presence of superfluidity, and scales as the number of 
atoms in the target. This provides an easily identifiable 
signature of the many body phase. 

Matter waves offer an effective method for non- 
destructively probing the many body phase in an optical 




FIG. 2: Cross sections for a probe with initial energy 
(Hk ) 2 /{2m) = E r , a lattice of depth V a = 15E r , with 4 
atoms and 4 sites, and J = 0.002. The values of U/J shown 
are: 0, 4, 8, and 16. The superfluid analytic result (light gray) 
for the inelastic background scales as N. The corresponding 
background for the Mott insulator (dark gray) is strongly sup- 
pressed. The inset shows the numerical cross section at the 
vertical dashed line as a function of U/J. 



lattice. More broadly, the scattering of matter waves de- 
pends strongly on properties such as the distribution of 
atoms in a partially filled lattice, and on the long range 
correlations and other manifestly quantum mechanical 
effects found in this novel state of matter. This makes 
this technique well-suited for examining a range of phe- 
nomena, including the effect of inhomogeneity of bosons 
and fermions and disorder in the lattice, in addition to 
probing the many body quantum phase transition. 
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